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Abstract 

We study the two-baryon system with two units of charm looking for the possible existence of 
a loosely bound state resembling the H dibaryon. We make use of a chiral constituent quark 
model tuned in the description of the baryon and meson spectra as well as the NN interaction. 
The presence of the heavy quarks makes the interaction simpler than in light baryon systems. 
We analyze possible quark-Pauli effects that would be present in spin-isospin saturated channels. 
Our results point to the non-existence of a charmed H-like dibaryon, although it may appear as a 
resonance above the AcAc threshold. 

PACS numbers: 12.39.Pn,14.20.-c,12.40.Yx 
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I. INTRODUCTION 


Hadronic molecules have recently become a hot topic on the study of meson spectroscopy 
as well-suited candidates for the structure of the so-called XYZ states [T] . This explanation 
seems compulsory in the case of the charged charmonium-like states hrst reported by Belle |2] 
and charged bottomonium-like states reported later on also by Belle [3]. In the two-baryon 
system there is a well-established molecule composed of two light baryons, the deuteron. 
Another well-known candidate is the H dibaryon suggested by Jaffe |1], a bound state with 
quantum numbers of the AA system, i.e., strangeness —2 and (T)J^ = (O)O'''. 

These two scenarios have recently triggered several studies about the possible existence of 
molecules made of heavy baryons [5H9]. The main motivation originates from the reduction 
of the kinetic energy due to large reduced mass as compared to systems made of light 
baryons. However, such molecular states that have been intriguing objects of investigations 
and speculations for many years, are usually the concatenation of several effects and not 
just a fairly attractive interaction. The coupling between close channels or the contribution 
of non-central forces used to play a key role for their existence. Some of these contributions 
may be reinforced by the presence of heavy quarks while others will become weaker. 

The understanding of the hadron-hadron interaction is an important topic nowadays. To 
encourage new experiments seeking for evidence of theoretical predictions, it is essential to 
make a detailed theoretical investigation of the possible existence of bound states, despite 
some uncertainty in contemporary interaction models ra- It is the purpose of this work 
to analyze the interaction of two-baryons with two units of charm and its application to 
study the possible existence of a hadronic molecule with the quantum numbers of the A^A^ 
system, {T)J^ = (O)O’''. When tackling this problem, one has to manage with an important 
difficulty, namely the complete lack of experimental data. Thus, the generalization of models 
describing the two-hadron interaction in the light flavor sector could offer insight about 
the unknown interaction of hadrons with heavy flavors. Following these ideas, we will 
make use of a chiral constituent quark model (CCQM) tuned on the description of the 
NN interaction m as well as the meson [T2] and baryon [THl [T3] spectrum in all flavor 
sectors, to obtain parameter-free predictions that may be testable in future experiments. 
Such a project was already undertaken for the interaction between charmed mesons with 
reasonable predictions na, what encourages us in the present challenge. Let us note that 
the study of the interaction between charmed baryons has become an interesting subject in 
several contexts [I6] and it may shed light on the possible existence of exotic nuclei with 
heavy flavors na. 

The paper is organized as follows. We will use Sec. [IT] for describing all technical details 
of our calculation. In particular. Sec. IIA presents the description of the two-baryon sys¬ 
tem quark-model wave function, analyzing its short-range behavior looking for quark-Pauli 
effects 


Sec. 


HB| briefly reviews the interacting potential, and section |II C| deals with the 

we 


solution of the two-body problem by means of the Fredholm determinant. In Sec. HI 


present our results. We will discuss the baryon-baryon interactions emphasizing those as¬ 
pects that may produce different results from purely hadronic theories. We will analyze the 
character of the interaction in the different isospin-spin channels, looking for the attractive 
ones that may lodge resonances. We will also compare with existing results in the literature. 
Finally, in Sec. |IV| we summarize our main conclusions. 
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II. THE TWO-BARYON SYSTEM 


A. The two-baryon wave function 


We describe the baryon-baryon system by means of a constituent quark cluster model, 
i.e., baryons are described as clusters of three constituent quarks. Assuming a two-center 
shell model the wave function of an arbitrary baryon-baryon system can be written as: 




LST 

B1B2 


(R) 



+ (-It 



( 1 ) 


where A is the antisymmetrization operator accounting for the possible existence of identical 
quarks inside the hadrons. The symmetry factor / satishes L + Si + S 2 — S+Ti+T 2 — T + f = 
odd m- For non-identical baryons indicates the symmetry associated to a given set of values 
LST. The non-possible symmetries correspond to forbidden states. For identical baryons, 
Bi = B 2 , f has to be even in order to have a non-vanishing wave function, recovering the 
well-known selection rule L -|- S' -|- T = odd. In the case we are interested in, two baryons 
with a charmed quark, the antisymmetrization operator comes given by 


A 




i=l j=4 


( 2 ) 


where B^^^ exchanges a pair of identical quarks i and j and P exchanges identical baryons. 
There appear two different contributions coming either from the exchange of light quarks 
(i = 1,2 and j = 4,5) or the exchange of the charm quarks (i = 3 and j = 6). If we 
assume gaussian Os wave functions for the quarks inside the hadrons, the normalization of 
the two-baryon wave function (i?) of Eq. (1) can be expressed as. 


- Ci(s,r) - c,(s,T)Ar,UR). 


(3) 





FIG. 1: Different diagrams contributing to the two-baryon normalization kernel as indicated in 
Eq. Q. The vertical thin solid lines represent a light quark, u or d, while the vertical thick solid 
lines represent the charm quarks. 
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A/"^(-R), Af^i{R), and A/"ex 2 (-^) stand for the direct and exchange radial normalizations de¬ 
picted in Fig. whose explicit expressions are 


Km 

Km 

KAR) 


dvr exp 
dvr exp 
dTrexp 



h+i/2 

*L-|-l/2 

h+l/2 



(4) 


where, for the sake of generality, we have assumed different gaussian parameters for the 
wave function of the light quarks (6) and the heavy quark [he). In the limit where the two 
hadrons overlap (i? —)■ 0), the Pauli principle may impose antisymmetry requirements not 
present in a hadronic description. Such effects, if any, will be prominent for L = 0. Using 
the asymptotic form of the Bessel functions, fL+ 1 / 2 , we obtain. 




rWa 2 


di —t dvr <1 — I To 4" I 2 

R^o 1 8 \b^ hi 


^ 6 V 8 [bK b2 


+ 


urn 


R^ f A 2 


exl -^ dvr < 1-^ TX + IT 

R^o 1 8 V 52 bl 


1 ^fRK‘^ 

^ 6 ( 8 [bl 


+ ... 


(5) 


Afm 


R^ f A 2 


-2 i ^ 8 I 62 + 52 


^ 6 I 8 bl 


+ 


Finally, the normalization kernel of Eq. (|^ can be written for the S wave in the overlapping 
region as, 

IS - f (i +1)} (ii - ^>1 + ■■■>■ («> 

where C{S,T) = |Ci(S',T) -|- 4C2(*S',T) for systems with two charmed baryons. For the 
particular case of the NEcc system C{S,T) = ^Ci{S,Ty. The values of C{S,T) are given 
in Table Thus, the closer the value of C{S,T) to 1/3 the larger the suppression of the 
normalization of the wave function at short distances, generating Pauli repulsion. It is 
the channel AcAc (T, J) = (0,0), with C{S,T) = 2/9^ where the norm kernel gets smaller 


TABLE I: C{S,T) spin-isospin coefficients, as defined in the text, for L = 0 partial waves. 


BiBj 

T = 0 

AcAc N 

T = 1 

R^cc AgEc EgEc 

T = 2 

EcEc 

II II 

0 

2/9 -1/9 0 

- -2/27 - 

5/27 -1/18 - 

16/81 1/54 -5/81 

-1/9 


^ Note that in the iVSec case the antisymmetrization operator of Eq. ([^ becomes much more simpler 
yf = (1 — 3Py^^), and only the first exchange diagram of Fig. will contribute with the two charm 
quarks on the same baryon. 
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FIG. 2: Normalization kernel as defined in Eq. Q for L = 0 and four different channels. 

at short distances. One would only find Pauli blocking [TH] in excited states like 
(T,J) = (2,3), where C{S,T) = 1/3, due to lacking degrees of freedom to accommodate 
the light quarks present on this conhguration, four u quarks with spin up. However, this 
partial wave could only exist for L = odd, and then the Pauli blocking may get masked by 
the centrifugal barrier. 

We show in Fig. |^the normalization kernel given by Eq. (|^ for L = 0 and four different 
channels: AcAc with (T, J) = (0,0), ScEc with (T, J) = (0,0) and (2,0), and NEcc with 
iT,J) = (1,0). In the hrst case C{S,T) is positive and close to 1/3, what gives a small 
normalization kernel. In the last case C{S,T) is also positive but smaller, giving rise to 
a slightly larger normalization kernel. In the other two cases C{S,T) is zero or negative, 
showing a large norm kernel at short distances and therefore one does not expect any Pauli 
effect at all. 


B. The two-body interactions 

The interactions involved in the study of the two-baryon system are obtained from a chiral 
constituent quark model m- This model was proposed in the early 90’s in an attempt to 
obtain a simultaneous description of the nucleon-nucleon interaction and the light baryon 
spectra. It was later on generalized to all flavor sectors HZ], In this model hadrons are 
described as clusters of three interacting massive (constituent) quarks, the mass coming 
from the spontaneous breaking of the original SU{2)l G) SU{2)ji chiral symmetry of the 
QCD Lagrangian. QCD perturbative effects are taken into account through the one-gluon- 
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exchange (OGE) potential [19]. It reads, 


VoGE{rij) 


a.y ?,[ 1 1 M I 1 I 2ai-a, \ 3S„ ] ^ 

4 * I rij A\2m‘f 2 m^ Sruimj J Vij Arriimjrfj \ 


(7) 


where are the SU{3) color matrices, ro = fo//i is a flavor-dependent regularization scaling 
with the reduced mass of the interacting pair, and as is the scale-dependent strong coupling 
constant given by 1121 . 

” In + ^)hl]' ® 

where oq = 2.118, no = 36.976 MeV and 70 = 0.113 fm“h This equation gives rise to 
as ~ 0.54 for the light-quark sector, as ~ 0.43 for uc pairs, and as ~ 0.29 for cc pairs. 

Non-perturbative effects are due to the spontaneous breaking of the original chiral sym¬ 
metry at some momentum scale. In this domain of momenta, light quarks interact through 
Goldstone boson exchange potentials. 


7x(Gj) — Gose(Gj) + GoPE(h'p), 


(9) 


where 


hbsE(f^p) 

VbpE(f^ ij) 


_9k 

dvr — mi 

qIe ml 




A2 


Y{ma Tij) - —Y{A Tij) 


dvr 12mimj A^ — 


rrir. 


A3 

Yirn^nj) - —Y{Anj) 


m-i 


CTi ’ (Jj 


A3 

H{m^rij) - —H{Arij) 


m-t 


Sij / (t^ * Td) . 


( 10 ) 


is the chiral coupling constant, Y{x) is the standard Yukawa function defined by 
Y(x) = jx, Sij = 3 ((jj ■ ■ fij) — <7* ■ aj is the quark tensor operator, and H{x) = 

(1 -I- 3/a; -I- 3/x^) Y{x). 

Finally, any model imitating QGD should incorporate confinement. Being a basic term 
from the spectroscopic point of view it is negligible for the hadron-hadron interaction. Lat¬ 
tice calculations suggest a screening effect on the potential when increasing the interquark 
distance 1201, 

VcoN(n,) = {-a, (1 - • A',). (11) 

Once perturbative (one-gluon exchange) and nonperturbative (confinement and chiral sym¬ 
metry breaking) aspects of QGD have been considered, one ends up with a quark-quark 
interaction of the form 


f hQj = nn] => VcoN(rij) + Voge(Ap) + 
I [QiQj = cn/cc] ^ V'con(Ap) V'oge(Ap) 


where n stands for the light quarks u and d. Notice that for the particular case of heavy 
quarks (c or b) chiral symmetry is explicitly broken and therefore boson exchanges do not 
contribute. The parameters of the model are the same that have been used for the study of 
the ND system |2l] and for completeness are quoted in Table [ll} The model guarantees a 
nice description of the light |T3| and charmed jTT] baryon spectra. 
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In order to derive the BnBm —>■ B^Bi interaction from the basic qq interaction dehned 
above, we use a Born-Oppenheimer approximation. Explicitly, the potential is calculated as 
follows, 

yB„B^{LST)^B^Bi{L> S'T){R) = ( 13 ) 

where 




iuL'S'T 
^BkBi 


m I E?<,=i 1^. 




vpLST 
r> R 

-On -O 777 


{R) 


^^bS^{R) 


S'T 

^BkBi 


(R) 


\hLST 

R R 


(R) 


nUiR) 


(14) 


In the last expression the quark coordinates are integrated out keeping R hxed, the result¬ 
ing interaction being a function of the two-baryon relative dis tance . The wave function 
^B^Bm(R) two-baryon system has been discussed in Sec. IIA 


C. Integral equations for the two-body systems 

To study the possible existence of two-baryon molecular states with two units of charm: 
AcAc, NScc, AcSc, or ScSc, we have solved the Lippmann-Schwinger equation for negative 
energies looking at the Fredholm determinant Dp{E) at zero energy [22]. If there are no 
interactions then Dp{R) = 1, if the system is attractive then Dp{R) < 1, and if a bound 
state exists then /li?(0) < 0. This method permits us to obtain robust predictions even for 
zero-energy bound states, and gives information about attractive channels that may lodge 
a resonance in similar systems [T^. We consider a two-baryon system BiBj in a relative S 
state interacting through a potential V that contains a tensor force. Then, in general, there 
is a coupling to the B^Bj D wave. Moreover, the two-baryon system can couple to other 
two-baryon states. We show in Table [TTTjthe two-baryon coupled channels in the isospin-spin 
(T, J) basis. 

Thus, if we denote the different two-baryon systems as channel A,, the Lippmann- 
Schwinger equation for the baryon-baryon scattering becomes 

7=Ai,A2,-" -^7=0,2 

X , a,/? = Ai, As, • • • , (15) 

where t is the two-body scattering amplitude, T, J, and E are the isospin, total angular 
momentum and energy of the system, iaSa, and are the initial, intermediate, and 


pldp-fV^ 


^q;Sq;)'^7<57 


7 ^ OL'y]T J 


{Pa,P'y) 


TABLE II: Quark-model parameters. 


mu,d{MeY) 

313 

5c\/(4vr) 

0.54 

mc(MeV) 

1752 

mo-(fm“^) 

3.42 

6(fm) 

0.518 

m,r(fm~^) 

0.70 

6c(fm) 

0.6 

A(fm“^) 

4.2 

ro (MeV fm) 

28.170 

Oc (MeV) 

230 

He (fm-Q 

0.70 
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final orbital angular momentum and spin, respectively, and is the relative momentum of 
the two-body system 7. The propagators G-y{E;p^) are given by 


p^') 


2 /i- 


^-pi^ + ie' 


with 


E= ^ 


2/^7 ’ 


(16) 


(17) 


where is the reduced mass of the two-body system 7. For bound-state problems E < Q 
so that the singularity of the propagator is never touched and we can forget the ie in the 
denominator. If we make the change of variables 


p^ = d 


1 + x-^ 

1 — x^’’ 


(18) 


where d is a scale parameter, and the same for pa and p^, we can write Eq. (15) as 

2 


^ — A. ... 0 —n 0^—1 \ 


X ^7) G,(E;p,) E). 


7=Ai,A2,-" -^7=0,2 ' 

S-y 




2d 

(1 - 

(19) 


We solve this equation by replacing the integral from —1 to 1 by a Gauss-Legendre quadra¬ 
ture which results in the set of linear equations 


N 


E \ ^ \ ^ j,jn£aSa,m£jSj 

7=Ai,j42,--- £ 7 = 0,2 m=l 


{E) (^-1 xu,E) = {Xr,, X ,), (20) 


with 


m: 


niaSa 


a7;T J 


(-^) SaS^ '^md 


-C 7 J/~v o 7 

/■£q; Sck r'^7 ^7 


l + Xr^ 


Xr 


2d 


(1 


Xr 


X y77‘-'(x„,a:,„)G7£;;p,„), (21) 

and where Wm and Xm are the weights and abscissas of the Gauss-Legendre quadrature while 


p^^ is obtained by putting x^ = Xm in Eq. (18). If a bound state exists at an energy Eb-, 

irTL^fy S jTTl-^xy S-y 

^a7;TJ 


the determinant of the matrix {Eb) vanishes, i.e., \Ma^.^Tj{EB) \ = 0. 


TABLE III: S and D wave two-baryon channels contributing to the different isospin-spin (T, J) 
states. See text for details. 
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FIG. 3: Different two-body potentials contributing to the {T)J^ = (0)0'*' channel. 


III. RESULTS AND DISCUSSION 


We will first discuss the interactions derived with the CCQM, centering our attention in 
the most interesting channel, the flavor singlet channel (T)J^ = (O)O^ with the quantum 
numbers of the AcAc state. This channel might lodge a charmed AT-like dibaryon. We show 
in Fig. I^the diagonal and transition central potentials contributing to the (T)J^ = (O)O’'' 
state. It is important to note that the AcAc system is decoupled from the closest two-baryon 
threshold, the state, that in the case of the strange H dibaryon becomes relevant for 
its possible bound or resonant character [2S]- The binding of the (T)J^ = (0)0+ state would 
then require a stronger attraction in the diagonal channels or a stronger coupling to the 
heavier ScSc state, that as we will discuss below is not fulfilled. 

In Fig. 1^ we have separated the contributions of the different terms in Eq. (12) to the two 
diagonal interactions. As can be seen, the AcAc potential is the most repulsive one. It be¬ 
comes repulsive at short-range partially due to the reduction of the normalization kernel (see 
Fig. §. The OGE and OPE can only give contributions through quark-exchange diagrams 
due to the color-spin-isospin structure of the antisymmetry operator [II1I21]. They gener¬ 
ate short-range repulsion that it is compensated at intermediate distances by the attraction 
coming from the scalar exchange, with a longer range. Thus, the total potential becomes 
slightly attractive at intermediate distances but repulsive at short range. In the EcEc in¬ 
teraction, the presence of a direct (without simultaneous quark-exchange) contribution of 
the OPE and the opposite sign of most part of the exchange diagrams, generates an overall 
attractive potential. This is rather similar to the situation in the strange sector but with the 
absence of the one-kaon exchange potential, what gives rise to a less attractive interaction. 
Regarding the character of the interaction, similar results were obtained in Ref. [9] within 
the quark delocalization color screening model (QDCSM). It is important to note at this 
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point the difference with hadronic potential models as those of Refs. EE]. As can be seen 
in Fig. 3(a) of Ref. [5j, the (T)J^ = (O)O’'' AcAc potential is attractive due to the absence 
of quark-exchange contributions and the dominance of the attraction of the scalar exchange 
potential. In spite of being attractive, the central potential alone is not enough to generate 
a bound state. In Ref. [B] they only consider the hadronic one-pion exchange and then, the 
AcAc interaction is zero. Thus all possible attraction comes generated by the coupling to 
larger mass channels. 

As mentioned above, when comparing with the similar problem in the strange sector an 
important difference arises, the absence of the AcAc O NScc coupling. As a consequence 
the mass difference between the two coupled channels in the {T)J^ = (O)O^ partial wave, 
AcAc and ScSc, is much larger than in the strange sector, making the coupled channel effect 
less important. Let us note that in the strange sector M(A^S) — M(AA) = 25 MeV and 
M(SS) — M(AA) = 154 MeV, this is why the NS channel plays a relevant role for the 
AA system |23], as well as why the NS state is relevant for the NA system [25]. In the 
charmed sector the closest channel coupled to AcAc in the (T)J^ = (O)O’'' state is ScSc, 338 
MeV above. This energy difference is similar to the NA — NN mass difference, the coupled 
channel effect being still important although it may not proceed through the central terms 
due to angular momentum selection rules [26|. Heavier channels play a minor role, as it 
occurs with the A A channel, 584 MeV above the NN threshold [26]. In the present case 
the coupling to the closest channel ScSc proceeds through the central potential and one 
does not expect higher channels, as S*S* 468 MeV above the AcAc threshold, to play a 
relevant role in quark-model descriptions as shown in the QDCSM model of Ref. [B]. The 




FIG. 4: Left panel: Contribution of the different terms of the interaction to the {T)J^ = (0)0"^ 
AcAc potential. ’OGE’ stands for the one-gluon exchange, ’OPE’ for the one-pion exchange, ’OSE’ 
denotes the one-sigma exchange and ’Tot’ represents the total potential. Right panel: Same as the 
left panel for the {T)J^ = (O)O''' ScNc potential. 
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FIG. 5: Left panel: Fredholm determinant of the {T)J^ = (O)O^ channel. The dashed line only 
considers the AcAc state, whereas the solid line includes the coupling to the ScSc state. Right 
panel: Fredholm determinant of the (T) = (0)0'*' coupled channel state for different values of be 

in fm. See text for details. 


situation seems to be a bit different in hadronic models where the non-central potentials 
are not regularized by the quark-model wave function Let us hnally note that the 

coupling between the AcAc and SeSe channels comes mainly given by quark-exchange effects 
and the direct one-pion exchange potential. Thus, it becomes a little bit stronger than in 
hadronic theories based on the one-pion exchange potential [6]. The resulting interaction is 
rather similar to that in the strange sector, as can be seen by comparing with Fig. 1(b) of 
Ref. |2Z], the main difference coming from the behavior of the normalization kernel at short 
distances. 

With these ideas in mind and following the procedure described in Sec. IIC we have 
performed a full coupled channel calculation of the (T)J^ = (O)O^ state. The results are 
shown in Fig. In the left panel we show by the dashed line the Fredholm determinant 
of the AcAc channel alone. The Fredholm determinant is large, indicating a barely small 
attractive interaction. When the coupling to the heavier SeSe channel is included, solid line 
in Fig. the system gains attraction, but it is not enough as to get a bound state. The main 
uncertainty when determining the baryon-baryon interaction in quark models with charmed 
baryons would be the harmonic oscillator parameter of the charm quark. We have explored 
the results for different values of be- The results are shown in the right panel of Fig. |^and 
as can be seen in no case the (T)J^ = (0)0+ state would become bound. Note that in 
Ref. [28] it was argued that the smaller values of be are preferred to get consistency with 
calculations based on inhnite expansions, as hyperspherical harmonic expansions [29] , where 
the quark wave function is not postulated. This also agrees with simple harmonic oscillator 
relations be = bn W—. The smaller values of be give rise to the less attractive results. For the 
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larger values of be, if a loosely bound state could be generated, the electromagnetic repulsion 
arising in the {T)J^ = (O)O’'' channel due to the electric charge of the A+ might dismantle 
the bound state. 

Thus, without the strong transition potentials reported in the QDCSM model of Ref. [S] or 
the strong tensor couplings occurring in the hadronic one-pion exchange models of Refs. m, 
it seems difficult to get a bound state in this system. We have recently illustrated within the 
quark model isni how the coupled channel effect between channels with an almost negligible 
interaction in the lower mass channel works for generating bound states. Although for the 
four-quark problem, in this reference it is demonstrated (see Fig. 2 of Ref. |30]) how when 
the thresholds mass difference increases, the effect of the coupled channel diminishes, which 
is an unavoidable consequence of having the same hamiltonian to describe the hadron masses 
and the hadron-hadron interactions. _ 

We have also analyzed the other {T)J^ channels shown in Table III with similar conclu¬ 
sions, the weak interaction in the charm sector and the absence of channel coupling between 
close mass channels works against the possibility of having dibaryons with two units of 
charm. For the sake of completeness we have calculated the Fredholm determinant for all 
cases and it is shown in Fig. 

One should hnally note that the problem of double heavy dibaryons has also been ap¬ 
proached in the literature by means of six-quark calculations. The group of Grenoble [31] 
addressed this problem within a pure chromomagnetic interaction obtaining several candi¬ 
dates to be bound. There also interesting results based on relativistic six-quark equations 
constructed in the framework of the dispersion relation technique [32] with a rich spec¬ 
troscopy of double charmed and beauty heavy dibaryons. Future experimental results will 
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-10 -5 0 


E (MeV) 



FIG. 6: Left panel: Fredholm determinant of the {T)J^ channels where the N'Ecc state is the 
lowest threshold. Right panel: Fredholm determinant of the {T)J^ channel where the SeSe state 
is the lowest threshold. 
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help to scrutinize among the different models, and in this way to improve our phenomenolog¬ 
ical understanding of QCD in the highly non-perturbative low-energy regime. This challenge 
could only be achieved by means of a cooperative experimental and theoretical effort. 


IV. SUMMARY 

In short summary, we have studied the baryon-baryon interaction with two units of charm 
making use of a chiral constituent quark model tuned in the description of the baryon and 
meson spectra as well as the NN interaction. Several effects conspire against the existence 
of a loosely bound state resembling the H dibaryon. First, the interaction is weaker than 
in the strange sector. Second, there is no coupling between close mass thresholds, like 
AcAc y-)- NScc, the closest threshold being more than 300 MeV above. Finally, the existence 
of a weak attraction may be killed by the electromagnetic repulsion absent in the strange 
sector. Thus, our results point to the nonexistence of low-energy dibaryons with two units 
of charm and in particular, to the nonexistence of a stable charmed H-like dibaryon. Given 
that the interaction in the {T)J^ = (O)O’'' is attractive, this state may appear as a resonance 
above the A^A^ threshold. 

Weakly bound states are usually very sensitive to potential details and therefore the¬ 
oretical investigations with different phenomenological models are highly desirable. The 
existence of these states could be scrutinized in the future at the LHC, J-PARC and RHIC 
providing a great opportunity for extending our knowledge to some unreached part in our 
matter world. 
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